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Sucrose and trehalose-6-phosphate (T6P) are central compounds in the regulation and 
orchestration of whole plant metabolism, growth, development, and flowering. To evaluate 
their highly complex and regulatory interaction with the two conserved sugar and energy 
sensors Snf1-related protein kinase 1 (SnRKI), an AMPK-related protein kinase, and 
hexokinase (Hxk), we developed a kinetic model which demonstrates the subtle metabolic 
control of sugar homeostasis in a wide range of concentrations without the need for 
changes in gene expression or protein concentrations. Our model approach is based on a 
comprehensive set of published metabolite concentrations under various conditions and 
coupled enzyme kinetics accounting for the role of SnRKI and Hxk in the sugar and energy 
homeostasis. This allowed us to investigate interactions between sugar phosphates, such 
asTBR which are metabolic inhibitors of SnRKI and Hxk, and sucrose synthesis during the 
transition from carbon deficiency to availability. Model simulations and sensitivity analyses 
indicated that slight changes in SnRKI activity induced by allosteric effectors may be 
sufficient to explain a dramatic readjustment of metabolic homeostasis. This may comprise 
up to 10-fold changes in metabolite concentrations. Further, the Hxk/T6P/SnRK1 interaction 
implemented in the model supports the interpretation of phenotypic and transcriptomic 
changes observed in Hxk overexpressing plants. Finally, our approach presents a theoretical 
framework to kinetically link metabolic networks to underlying regulatory instances. 
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INTRODUCTION 

Plant carbohydrates are the primary products of photosynthe- 
sis and are thus central to regulation of metabolism, devel- 
opment and growth. Carbohydrates are synthesized as triose 
phosphates in chloroplasts. The carboxylation of ribulose-1,5- 
bisphosphate generates an instable C6 compound that decomposes 
into two molecules of glycerate- 3 -phosphate. Under consump- 
tion of NADPH and ATP being generated by the photosynthetic 
electron transport chain, glycerate- 3 -phosphate is reduced to give 
triose phosphates, which can either be exported to the cytosol 
or be metabolized in the chloroplast as a substrate for starch 
synthesis (Heldt etal., 2011). Cytosolic triose phosphates are sub- 
strate for sucrose synthesis for which fructose- 1,6-bisphosphatase 
(cFBPase) and sucrose phosphate synthase (SPS) have been iden- 
tified as rate limiting (Stitt etal, 1983; Strand etal., 2000). 
Sucrose may then be transported to sink organs or it is cleaved 
by invertase to glucose and fructose (Sturm, 1999). Free hex- 
oses are phosphorylated in a reaction catalyzed by hexokinase 
(Hxk). The hexose phosphates are substrates for sucrose synthe- 
sis, glycolysis or the synthesis of organic and amino acids (Talts 
etal, 2004; Claeyssen and Rivoal, 2007; Ainsworth and Bush, 
2011). 

Beside their metabolic interaction via enzymatic interconver- 
sions, plant sugars act as signaling molecules. Analyses of the 



Arabidopsis glucose insensitive 2 mutant {gin2), which lacks the 
catalytic activity of a specific hexokinase (HXKl), have identi- 
fied AtHXKl as a central component in plant sugar sensing and 
signaling (Moore etal., 2003). Hxk is present in various cel- 
lular compartments, also comprising the nucleus (Yanagisawa 
etal, 2003). HXKl was suggested to form a glucose signal- 
ing complex core together with the vacuolar H+-ATPase Bl 
and the 19S regulatory particle of proteasome sub unit, RPT5B 
(Cho etal., 2006). In context of its metabolic activity, mito- 
chondrial Hxk together with mitochondrial and/or cytosolic 
invertase activity was discussed to contribute to the home- 
ostasis of reactive oxygen species (ROS; Camacho-Pereira etal., 
2009; Xiang etal., 2011). Finally, due to its connection and 
interaction with hormone signaling pathways (Rolland etal, 
2006; Hanson and Smeekens, 2009), it becomes obvious that 
Hxk plays various complex roles in metabolism of the whole 
plant. 

In addition to Hxk, there are several other central regulatory 
instances linking sugar metabolism to signaling networks and 
energy metabolism. One prominent example is the Snfl -Related 
Protein Kinase-1 (SnRKI; AKINlO/11) which has a catalytic 
domain similar to that of Snfl (Sucrose non-fermenting- 1 ) of yeast 
and AMPK (AMP-activated protein kinase) of animals (Hardie, 
2007). It has been shown that SnRKI is involved in the regulation 
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of a broad array of target genes which orchestrate the response 
of metaboHsm to starvation (Baena-Gonzalez etal., 2007). Cen- 
tral plant proteins identified as direct substrates for SnRKl are 
3-hydroxy-2-methylglutaryl-coenzyme A (HMG-CoA) reductase, 
nitrate reductase (NR), 6-phosphofructo-2-kinase/fructose-2,6- 
bisphosphatase (F2KP), TPS and SPS [for an overview see (Halford 
etal, 2003) and (Halford and Hey, 2009)]. SnRKl phosphoryla- 
tion results in inactivation of these enzymes, yet inactivation of 
NR, F2KP and TPS also requires the binding of a 14-3-3 protein 
(Bachmann etal, 1996; Moorhead etal, 1996; Halford and Hey, 
2009). 

Energy deprivation induces a comprehensive SnRKl -mediated 
transcriptional and metabolic reprogramming. Hence, it is not 
surprising that SnRKl activity is regulated on various levels of 
molecular organization. Recently, (Rodrigues etal., 2013) showed 
that two clade A type 2C protein phosphatases (PP2Cs) dephos- 
phorylated and inactivate SnRKl. As PP2Cs are repressors of the 
abscisic acid (ABA) pathway this indicates a way how cellular 
metabolism is connected to plant hormonal interactions affect- 
ing the whole plant metabolism, development and growth. An 
activation of SnRKl in Arabidopsis has been shown to be due to 
two upstream kinases, SnRKl -activating kinases (SnAK) 1 and 2 
(Crozet etal, 2010). 

Beyond the regulation of SnRKl by other kinases and phos- 
phatases, the inhibition of SnRKl by metabolites, and predomi- 
nantly by intermediates of the central carbohydrate metabolism, 
has been shown to have significant impact on its activity. The 
phosphorylated sugars trehalose-6-phosphate (T6P), glucose- 1- 
phosphate (GIP), glucose-6-phosphate (G6P) as well as ribose-5- 
phosphate (R5P) and ribulose-5-phosphate (Ru5P) were shown 
to reduce SnRKl activity significantly while fructose-6-phosphate 
(F6P) and uridine-5^-diphosphoglucose (UDPG) did not inhibit 
SnRKl (Nunes etal., 2013b). Particularly the T6P/SnRKl interac- 
tion has been focused in various recent studies. T6P was shown to 
play a crucial role in regulation of carbon utilization and growth 
in Arabidopsis thaliana (Schluepmann etal, 2003). Hence, like 
SnRKl, it represents a regulatory instance which is likely to exert 
major control on plant metabolism and development. In plants, 
T6P is synthesized from UDPG and G6P. This reaction step is cat- 
alyzed by trehalose phosphate synthase (TPS) and enzyme TPSl is 
discussed to account for most TPS activity in plants (Vandesteene 
etal., 2010). In a subsequent step, T6P is converted to trehalose 
(Tre) in a reaction catalyzed by trehalose phosphate phosphatase 
(TPP). Finally, trehalase catalyzes the cleavage of Tre in two Glc 
equivalents. Recently, TPSl was suggested to be required for timely 
initiation of flowering in Arabidopsis thaliana (Wahl etal., 2013). 
In another study, the reaction product of TPSl, T6P, was shown 
to respond significantly to external Sue supply when Sue was 
fed to seedlings (Lunn etal., 2006). Yet, Nunes and co-workers 
discussed T6P not to be a growth signal per se, but rather to 
prime gene expression through SnRKl for growth in response 
to Sue accumulation under sink-limited conditions (Nunes et al., 
2013a). 

Assembling this information about the interaction of the cen- 
tral carbohydrate metabolism and its interaction with SnRKl 
yields a complex picture of metabolic regulation. Even a sim- 
plification of this network which only accounts for the kinase 



interaction of SnRKl with SPS, which does not require the binding 
of a 14-3-3 protein for its inactivation, and including inhibitions 
exerted by metabolites on SnRKl but also on Hxk (Claeyssen and 
Rivoal, 2007) and invertase (Inv; Sturm, 1999), still yields a com- 
plex model of the (cytosolic) carbohydrate metabolism (Figure 1). 
Particularly under changing environmental conditions, for exam- 
ple the addition of sucrose to starved seedlings, it is hardly possible 
to intuitively draw a conclusive and comprehensive picture about 
regulatory consequences yielding the adjustment of a metabolic 
homeostasis. Focusing the question how SnRKl can function- 
ally be integrated in the regulation of the central carbohydrate 
metabolism in terms of enzyme kinetics, we present an approach 
of kinetic modeling which allows for the comprehensive analysis 
of the SnRKl /C-metabolism interaction during a repletion experi- 
ment of C-starved seedlings of Arabidopsis thaliana. Our approach 
accounts for the experimental findings of various studies, which 
have been published previously, and provides a representative and 
systematic overview of central metabolic consequences emanating 
from SnRKl and Hxk activity. 

RESULTS 

A KINETIC EQUATION FOR SnRKl ACTIVITY 

To estimate the impact of SnRKl on the central carbohydrate 
metabolism and, vice versa, to analyze how its activity is affected by 
dynamic changes of metabolic inhibitor concentrations, a kinetic 
equation was derived to simulate the rate vsnRKi by which SPS is 
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FIGURE 1 I Schematic overview of the cytosolic carbohydrate 
metabolism in Arabidopsis. Enzymatic interconversions are indicated in 
black arrows, enzymes tliat catalyze each step are written in blue. Red 
dashed lines indicate inhibitory interaction, green dashed lines indicate 
activation. SnRKl: Snfl-related protein kinase 1; PP: protein phosphatase; 
F16P: fructose-1,6-bisphosphate; F6P: fructose-6-phosphate; G6P: 
glucose-6-phosphate; GIP: glucose-1-phosphate; UDPG: 
uridine-5'-diphosphoglucose; S6P: sucrose-6-phosphate; T6P: 
trehalose-6-phosphate; Sue: sucrose; Tre: trehalose; Fru: fructose; Glc: 
glucose; SuCgxtemal- externally supplied sucrose; SucExport- sucrose 
exported to sinks; cFBPase: cytosolic fructose-1 ,6-bisphosphatase; PGI: 
phosphoglucoisomerase; PGM: phosphoglucomutase; UGPase: 
uridine-5'-diphosphoglucose pyrophosphorylase; SPS: sucrose phosphate 
synthase; TPS: trehalose phosphate synthase; SPP: sucrose phosphate 
phosphatase; TPP: trehalose phosphate phosphatase; Hxk: hexokinase; 
GIcK: glucokinase; FrcK: fructokinase; Inv: invertase; SuSy: sucrose 
synthase. 
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phosphorylated and inactivated (Eq. 1). 

VSnRKl — 

Vmax,SnRKl " [-EsPS.actl " (l + ^partial • :^~~") 
JgM,E.Ps... + [£sPS,ac.] .(1+^^ + ^^ + ^+ + 

ivi,_ct,pt,,act oro,ciciJ ^i^teP -^i.GlP J^i,G6V ^i,T6P --Ki.GlP -^i.TeP ' -Ki.GeP / 

Kinetic parameters, i.e., Vmax,SnRKl, ^partiab ^M,ESPS,act> J^i,T6P, 

^i,GiP and Ki,G6P were estimated and derived from literature (see 
Materials and Methods). £sps,act represents the concentration of 
activated dephosphorylated SPS peptide which was the substrate 
for SnRKl. As SPS is phosphorylated and inactivated by SnRKl at 
Ser-158 (Halford etal., 2003), the phosphorylation stoichiometry 
was 1. Simulation of vsnRKi under single and combined inhibitor 
concentrations ranging from 0 to 1 mM finally resulted in inhibitor 
kinetics (Figure 2) which were highly similar to the experimental 
findings reported by Nunes et al. (2013b). 

Increasing T6P induced a very fast decline of vsnRKi but was 
asymptotically stable at a value of about 28% of the Vmax,SnRKi- 
Increasing GIP levels resulted in a significant reduction of vsnRKi 
which was, however, not as pronounced as that of T6P. The 
increase of G6P showed least impact of these three metabolic 
inhibitors on vsnRKi which fully agrees with the experimental find- 
ings (Nunes etal., 2013b). Combinations of T6P, GIP, and G6P 
revealed the highest impact on vsnRKi by a T6P-G1P-G6P combi- 
nation resulting in a final activity of about 10% of the Vmax,SnRKi 
when all inhibitors had a concentration of 1 mM. This effect was 
the strongest, followed by T6P-G1P (-12% of Vmax,SnRKi), T6P- 

G6P (-18% of Vn,ax,SnRKl) and G1P-G6P (-28% of Vn,ax,SnRKl). 



These findings are also highly similar to the previously published 
experimental results. 

EMBEDDING OF SnRKl ACTIVITY IN THE CENTRAL CARBOHYDRATE 
METABOLISM 

A mathematical model, based on a system of ODEs, was pro- 
grammed according to the graphical model structure shown in 
Figure 1. A detailed description of the model and its design 
is provided in the section "Materials and Methods." The func- 
tional connection of SnRKl and SPS activity was arranged by 
two ODEs describing the time-dependent dynamics of concen- 
trations of phosphorylated (inactivated) and dephosphorylated 
(activated) SPS peptide, i.e., SPS protein with one phosphoryla- 
tion event per protein. Activated SPS peptides were substrate for 
SnRKl and simultaneously defined the maximum enzyme activity 
of SPS (vmax,sps) as described in Eq. 2. Model simulations aimed at 
the reproduction of a short-term sugar feeding experiment with 
C-starved Arabidopsis seedlings (Yadav etal, 2014). Hence, the 
simulation time equaled the time period of external supply of 
C-starved seedlings with 15 mM sucrose. The initial conditions 
were given by the metabolite levels of seedlings before the start 
of sucrose supply. Except for trehalose, activated and inactivated 
SPS peptide all initial concentrations were taken from the study 
of Yadav etal. (2014). The initial trehalose content was estimated 
to be in the range of nM/low |xM (Paul etal., 2008). The inacti- 
vated and activated SPS peptide amount was estimated from data 
reported by Lehmann etal. (2008): from a total amount of SPS 
protein, £sPS,total ^ 8-5 x 10~^ |xmol gFW~^, we assumed that 
C-starved seedlings compared to non-C-starved seedlings have a 




0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Inhibitor concentration [mlVI] 

FIGURE 2 I Simulated kinetic effects of metabolic inhibitors on SnRKl activity. Lines represent the rate of phosphorylation by SnRKl (^snRKl i-®-- "the 
solutions of equation (1), under various inhibitor concentrations and combinations. 
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higher activity of SnRKl, and hence the initial condition for acti- 
vated SPS, £sps,act) was determined to be significantly lower than 
for inactivated SPS, £sps,inact- Based on the comparison of SPS 
activity during day (high energy supply) and night (lower energy 
supply; Nagele etal, 2010, 2012), initial £sps,act of C-starved 
seedlings were estimated to be about 15-25% of £sPS,total result- 
ing in an amount of £sPS,act ranging 1.275 to 2.125 x 10~^ |xmol 
gFW-^. 

MODEL SIMULATIONS REVEAL A VERY FAST REGULATION OF SPS 
INDUCED BY THE DYNAMICS OF METABOLIC INHIBITORS OF SnRKl 

Simulating metabolite dynamics during an external sucrose sup- 
ply over 3 h resulted in metabolite dynamics (Figure 3) which 
were highly similar to previously published experimental find- 
ings (Yadav et al., 2014). All metabolite levels showed a significant 
increase during the first 30 min of external sucrose supply. Most of 
the metabolite pools increased following a saturating course, while 
particularly free hexoses, i.e., glucose and fructose, increased lin- 
early (Figures 3F,G) . To test whether our model and its underlying 
kinetic parameters and assumptions were also valid beyond a sim- 
ulation time of 3 h, we simulated an external 15 mM sucrose supply 
for 8 h which has been experimentally performed by Lunn et al. 
(2006). Although absolute metabolite levels differed between the 
two studies of Yadav et al. (2014) and Lunn et al. (2006) our model 
simulations resulted in the same qualitative metabolic changes 
as described by both of them: only hexose levels show a linear 
increase even after 8 h of external sucrose supply (Figure 3, small 
diagrams). These simulation results indicate that our proposed 
model successfully describes metabolite dynamics ranging from 
C-starved to sugar supplied seedlings. 

Simulation results further indicated the level of trehalose to be 
in a low range compared to other metabolite levels (Figure 3K). 
While we have no explicit experimental evidence for that simu- 
lation result, it is in agreement with the finding of Mliller et al. 
(2001) who could not determine trehalose in leaf tissue of Ara- 
hidopsis, indicating a very low absolute amount compared to other 
carbohydrates. 

A central output of our model simulations was the occurrence 
of a very fast and linear activation step of SPS which was due to 
an inactivation of SnRKl by the increasing levels of its metabolic 
inhibitors (Figure 4). Simulations revealed a time for activation 
of about 0.003 h or 11 s (Figure 4A), i.e., the first molecule 
of the external supply enters the intracellular reaction network 
and it takes about lis until £sps,act equals £sps,inact- The sec- 
ond activation step was slower and showed an asymptotic course 
(Figure 4B) finally resulting in about 60% activated SPS peptide 
(4.75 X 10"^ ixmolgFW"^ of 8.2 x 10"^ jxmolgFW-^). 

THE INTERACTION 0FT6P, SnRKl AND SUCROSE 

Current models for the role of T6P in plant growth and metabolic 
regulation discuss T6P to be closely related to sucrose avail- 
ability and to inhibit SnRKl which increases the expression of 
biosynthetic genes, e.g., Nunes etal. (2013a). In this context, 
we have analyzed the simulated T6P versus sucrose content and 
found, as already expected from the experimental data (Yadav 
etal, 2014), a corresponding picture (Figure 5A). Particularly 
in the sucrose concentration range 0.25-3.5 |xmol gFW~^, the 



T6P level increased exponentially while it became almost linear 
at higher sucrose concentrations (3.5-4.6 |xmol gFW~^). These 
higher sucrose levels emerged after at least 1.25 h of external 
sucrose apply (Figure 3E). At this time point, also the reaction 
rate of SnRKl became linear with respect to T6P concentrations 
(Figure 5B) indicating a tight connection of T6P and sucrose via 
SPS and SnRKl. 

To evaluate how the SnRKl kinetic is affected at different 
time points of external sucrose supply and under physiological 
conditions, the simulated concentrations of metabolic inhibitors 
were determined and vsnRKi was calculated (Figure 6). Due to 
the partial and non-competitive inhibition of SnRKl there is 
no change in its substrate affinity but only in the Vmax,SnRKi- 
These simulations indicate that a change from vsnRKi = 1-2 
before sucrose supply to vsnRKi = 1.0 |xmol Peptide h~^ gFW~^ 
after 3 h of sucrose supply can explain the observed metabolic 
changes by a reduced rate of SPS phosphorylation, i.e., inactiva- 
tion. 

STEADY-STATE SENSITIVITY ANALYSIS REVEALS A DIFFERENT 
BEHAVIOR OF GLUCO- AND FRUCTOKINASE DUE TO PERTURBATIONS 
IN EXOGENOUS SUGAR SUPPLY 

Due to the observed linear increase of free hexoses, which was 
also occurring after 8 h of simulation time, a steady-state of the 
system of ODEs, i.e., d/dt = 0 for all metabolite levels, could not 
be determined. There are several explanations for such a behavior 
and one of them is a missing element in the model structure, 
for example a subcellular compartment for metabolite storage 
like the vacuole. In a previous study, Kruger etal. (2007) have 
already demonstrated the impact of the vacuolar glucose pool 
on the calculation of steady-state fluxes in heterotrophic Ara- 
bidopsis tissues. In the vacuolar compartment, hexoses cannot 
be phosphorylated due to a missing hexokinase activity. At the 
same time, a different invertase isoenzyme with different biochem- 
ical properties is active (Sturm, 1999). Hence, we tested whether 
we could identify a metabolic steady- state after 3 h of external 
sucrose supply by variation of invertase-related parameters (vmax? 
Km, Ki). This approach succeeded and a metabolic steady- state 
could be identified with a reduced Vmax and an increased Km of 
invertase. 

In previous studies, hexokinase (Hxk) has been shown and 
discussed to be a central component in plant sugar signal- 
ing and perception [an overview is provided in (Rolland etal, 
2006)]. To analyse a possible interaction between the SnRKl - 
related and the Hxk-related regulatory influence on the central 
carbohydrate metabolism, a steady- state sensitivity analysis was 
performed focusing on changes which are induced by a per- 
turbation of the external sucrose supply. The regulatory impact 
of Hxk on the metabolic homeostasis was estimated by an in 
silico overexpression of Hxk activity, i.e., the Vmax of glucok- 
inase and fructokinase were increased threefold. In this way, a 
possible limitation of sugar metabolism by the hexose phosphory- 
lation rate would be significantly decreased and the susceptibility 
of other subsequent metabolic steps to a change in external 
sucrose supply could be identified. Ratios of sensitivities calcu- 
lated for the overexpressed and the wild type model are provided 
in Table 1. 
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FIGURE 3 I Simulation of metabolite dynamics induced by external sucrose supply. Results of simulations are provided for all metabolites (A-l) within a 
simulation period of 3 h (large diagrams) and 8 h (small diagrams). Except forTre (I), simulations of 3 h directly refer to experimental findings of Yadav etal. (2014). 



The comparison of steady- state sensitivities revealed an artic- 
ulate increase of sensitivities of vsps? UDPG, vprcK? as well as all 
compounds and reactions of the T6P/Tre metabolism due to the 
increased Vmax of gluco- and fructokinase. The strongest decrease 
in sensitivities was observed for vgIcK? ^inv and Glc. At this point, 



it is interesting to note that although both metabolic hexokinase- 
associated activities, gluco- and fructokinase, were simulated with 
a threefold increase in Vmax? the sensitivities of both reaction 
rates behaved differently - sensitivity of vprcK increased while the 
sensitivity of vgIcK decreased. 
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SimulationTime [h] SimulationTime [h] 

FIGURE 4 I Simulated time course of activated and inactivated SPS peptide. Black lines indicate activated/dephosphorylated SPS peptide, turquoise lines 
indicate inactivated/phosphorylated peptide. (A) shows the first 0.01 h after simulated external sucrose supply, (B) shows the whole simulation period of 3 h. 




Sucrose [|jnnol gFW""] T6P [|jmol gFW-^ ] 



FIGURE 5 I The SnRK1-T6P-Sucrose interaction. (A)TheT6P level first increases exponentially, then linear with the sucrose level. (B)The decrease of SnRKI 
activity with respect toT6P concentrations occurring during 3 h of external sucrose supply. 




FIGURE 6 I SnRKI kinetics under various physiological conditions. 

Simulations of SnRKI kinetics were performed solving equation (1) with the 
simulated levels of metabolic inhibitors as they are presented in Figure 3. 



Arrows on the top indicate the level of substrate, i.e., of dephosphorylated 
SPS peptide, at the time points 0 h (blue), 0.5 h (black), 1.25 h (red), and 3 h 
(green) after external sucrose supply (15 mM). 
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Table 1 | Changes in sensitivities to perturbation in external sucrose 
supply due to Hxk overexpression. 



Model states/reaction rates Sensitivity ratio Hxk ox vs WT 



vsps 


7.42 


UDPG 


5.83 


\/FrcK 


2.92 


T6P 


2.14 


Tre 


2.14 


vjps 


2.1 


Vjpp 


2.1 




2.1 


F6P 




Sue 




Fru 




^Suc. import 




ur D rcioc 




VpQ\ 




^Suc. Export 








^SPS.act 


0.93 


^SPS.inact 


0.9 


\/SnRK1 


n Q 
u.y 


Vpp 


0.9 


G6P 


0.18 


G1P 


0.18 


VPGM 


0.18 


\/UGPase 


0.18 




0.071 




0.07 


Glc 


0.061 



Ratios are sorted in a descendent order Vx: rate of tlie reaction catalyzed by 
enzyme x; all other metabolite and protein abbreviations refer to the legend in 
Figure 1. 

DISCUSSION 

Sensing and integrating environmental stimuli is a prerequisite for 
the successful adaptation of organisms to a changing environment. 
Due to their sessile lifestyle, particularly plants have developed 
various strategies to cope with environmental dynamics. Central 
and conserved mechanisms, such as the SNFl/AMPK/SnRKl pro- 
tein kinases and the sugar sensor hexokinase, have been shown 
to play a crucial role in regulation of metabolism, development 
and stress response (Rolland etal., 2006; Polge and Thomas, 
2007). In recent years, numerous comprehensive studies have 
focused on analyzing the role of SnRKl in the reprogramming 
of metabolism. SnRKl was identified to be significantly involved 
in linking stress, sugar, and developmental signals, thus playing a 
pivotal role in the global regulation of plant metabolism, energy 
balance and growth (Baena-Gonzalez etal, 2007). In addition, 
central metabolic signaling pathways comprising the sucrose and 
T6P metabolism have recently been discussed to have a major 



impact on metabolic and developmental homeostasis (Yadav et al, 
2014). Although these - and numerous other - studies provide 
compelling evidence for a highly efficient, interlinked and cen- 
tral regulatory network, its function in vivo and, particularly, 
its properties under physiological and dynamic conditions is still 
poorly understood. Beyond, it is still not clearly understood how 
SnRKl activity affects the metabolic reprogramming during stress 
conditions, for example C-starvation. While it is known that low 
levels of glucose, high levels of sucrose and darkness induce the 
activity of SnRKl (Rolland etal, 2006), a comprehensive charac- 
terization also accounting for the differential impact of T6P, GIP, 
and G6P on SnRKl activity (Nunes etal, 2013b) is difficult to 
obtain. This is mainly due to the complexity of the underlying 
metabolic and signaling networks comprising numerous compo- 
nents and a vast amount of regulatory interactions. During the 
last decade, mathematical kinetic modeling has emerged as an 
attractive approach to analyse such networks which produce out- 
puts that can hardly be assessed by intuition. Here, we applied 
such an approach to the sucrose/T6P metabolism by connect- 
ing SnRKl with SPS activity in a mathematical model based on 
enzyme kinetics. Although our model was solely based on allosteric 
and posttranslational modifications, and did neither include tran- 
scriptional regulation of genes encoding enzymes of carbohydrate 
metabolism nor activation/inactivation of SnRKl by SnAKs or 
PP2Cs, we could successfully simulate the experimental output 
from previous studies on C-starved Arabidopsis seedlings which 
were fed with sucrose (Lunn et al., 2006; Osuna et al., 2007; Yadav 
etal., 2014). While this model is far from being complete because 
it focuses on central interactions of proteins and metabolites and 
disregards various regulatory instances, it is remarkable that com- 
prehensive and significant metabolic changes can be accounted 
for. This indirectly supports the suggestion of Nunes et al. (2013a) 
who discussed T6P not to be a growth signal per se but to prime 
various processes, such as gene expression, through SnRKl. This 
indicates that, despite the compelling experimental evidence for 
multiple interlaced regulatory circuits involving the transcrip- 
tome, proteome and metabolome, the regulation of a metabolic 
homeostasis under certain physiological conditions may involve 
only a subset of the available regulatory instances. Finally, this 
again emphasizes the necessity of modeling approaches due to the 
diversity of explanations for metabolic regulation originating from 
the non-linear relationship between different levels of molecular 
organization. 

Our simulations show that a slight reduction ( <20%) of SnRKl 
activity, exerted by its metabolic effectors, is sufficient to explain a 
dramatic reprogramming of the central carbohydrate metabolism 
(Figure 6) . Additionally, sensitivity analysis provided evidence that 
the phosphorylation/dephosphorylation step of SPS was indis- 
pensable for a successful simulation (data not shown) because 
slight modification (1%) of SnRKl- and PP- related parameters 
had a dramatic effect on the model output. Finally, our simulations 
revealed that small variations in T6P concentration had a large 
impact on sucrose content (Figure 5A), particularly in the phys- 
iological range of 0.5-3 |xmol Sue gFW~^. This highly sensitive 
T6P-sucrose interaction sheds light on previous findings which 
showed T6P to be indispensable for carbohydrate utilization and, 
as a consequence, growth in Arabidopsis thaliana (Schluepmann 
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etal., 2003): sucrose represents the major transport sugar which 
is, at the same time, central to numerous regulatory processes 
(Koch, 2004) and stabilizes metabolism during environmental 
perturbation (Nagele etal, 2010). 

Further, it is interesting to note that above a T6P level of 0.3- 

0. 5 nmol gFW~^ and a sucrose level above 3.5 |xmol gFW~^ the 
relationship of both metabolites was almost linear (Figure 5A). 
In the same concentration range, also the SnRKl activity behaved 
linear with respect to T6P (Figure 5B), indicating an indirect regu- 
lation of sucrose by the T6P-SnRKl interaction. This might deliver 
the kinetic explanation for a previous finding, where increases of 
sucrose above this level were found to result in a proportionate 
increase both in T6P level and changes in expression of SnRKl 
marker genes (Nunes et al, 2013a). 

A major difference of allosteric regulation compared to tran- 
scriptional regulation is the time frame in which regulation finally 
takes place. Our calculations revealed that the initiation phase 
of SPS activation by reduction of its phosphorylation required 
10-15 s after external sucrose supply to C-starved seedlings 
which makes it a quite efficient regulatory system capable of a 
fast response to fluctuating environmental conditions. Hence, 
although T6P and the reaction rates involved in its metabolism 
represent only a fraction of those found in sucrose metabolism 
(Paul etal, 2008), we provide further evidence for its highly 
sensitive regulatory impact on sucrose metabolism. 

Addressing the question how two conserved central signal- 
ing mechanisms - T6P/SnrKl and Hxk - might be coordinated 
within the central carbohydrate metabolism, we performed an in 
silico overexpression of Hxk, i.e., GlcK and FrcK activity. Hex- 
okinase activity represents a rate limiting step in sucrose cycling, 

1. e., the cleavage and re-synthesis of sucrose. Previously we could 
show that it participates significantly in the stabilization of the 
central carbohydrate metabolism due to environmental perturba- 
tions (Nagele et al., 2010; Henkel et al., 201 1). Beyond its catalytic 
activity, HXKl was shown to be a glucose sensor integrating nutri- 
ent, light, and hormonal signaling (Moore etal., 2003). Hence, 
in silico overexpression of Hxk activities increased the maximum 
capacity of hexose phosphorylation, thus reducing its rate lim- 
itation. The finding that, as a consequence of increased Hxk 
activities, the whole T6P metabolism was affected more than twice 
as much as in the wild type model under external perturbations 
(Table 1) indicated a possible interdependency of Hxk- and T6P- 
dependent signaling pathways. Additionally, the strongest effect of 
environmental perturbation was observed for sucrose synthesis, 
catalyzed by SPS, which is directly related to the T6P sensitivi- 
ties by SnRKl. Yet, as the phosphorylation rate of SnRKl, vsnRKi? 
was only slightly affected in its sensitivity, the main perturba- 
tion of the T6P/SnRKl/SPS interaction might rather be related 
to the UDPG pool. Due to being a substrate for both SPS and 
TPS, changes in the UDPG pool as well as in its sensitivity 
and stability are also affecting the pathways of sucrose and T6P 
synthesis. 

Based on our findings we hypothesize that an increase of Hxk 
activity, e.g., by AtHXKl overexpression, might induce fluctua- 
tions in the T6P metabolism which can directly, and very fast, affect 
SnRKl activity and/or other signaling compounds which have not 
been analyzed here. In this context, previous studies reporting 



on AtHXKl overexpression in tomato and Arabidopsis found a 
decreased expression of photosynthetic genes and an inhibition of 
growth (Dai et al., 1999; Kelly et al., 2012). This physiological out- 
put is similar to the output being characteristic for SnRKl activity 
(Baena- Gonzalez etal., 2007), and hence it supports a model of 
Hxk/T6P/SnRKl -interdependent signaling. 

MATERIALS AND METHODS 
MODEL DRAFT 

A model of the central cytosolic carbohydrate metabolism in 
plants was conceptualised comprising the enzymatic intercon- 
version of the metabolite pools of fructose 6-phosphate (F6P), 
glucose 6-phosphate (G6P), glucose 1 -phosphate (GIP), uri- 
dine 5-diphosphoglucose (UDPG), trehalose 6-phosphate (T6P), 
trehalose (Tre), sucrose 6-phosphate (S6P), sucrose (Sue), glu- 
cose (Glc) and fructose (Fru). One model input was defined 
by the rate of interconversion from fructose 1,6-bisphophate 
(F16P) to F6P which is catalyzed by cytosolic fructose 1,6- 
bisphosphatase (cFBPase). cFBPase activity was shown to be 
rate limiting for sucrose biosynthesis (Strand etal., 2000), and 
hence subsequent reactions catalyzed by the enzymes phospho- 
glucoisomerase (PGI), phosphoglucomutase (PGM) and UDPG 
pyrophosphorylase (UGPase) were assumed to be not rate lim- 
iting (Meng etal., 2009). We are aware that this represents a 
strong simplification of the in vivo system, which, in addition 
to cFBPase, also comprises the enzymatic activity of inorganic 
pyrophosphate-dependent diphosphate-fructose-6-phosphate 1- 
phosphotransferase (PFP) and ATP-dependent phosphofructok- 
inase (PFK). Yet, as our approach focuses on the subsequent 
steps of the cFBPase/PFP/PFK-driven reaction, we omitted these 
reaction kinetics and the regulatory impact of the fructose- 
6-phosphate 2-kinase/fructose-2,6-bisphosphate 2-phosphatase 
reaction to reduce uncertainties in the biochemical model being 
caused by structure, parameters and kinetics (Schaber et al., 2009) . 
To simulate the effect of external sucrose on the central car- 
bohydrate metabolism, a second model input was defined by a 
Michaelis-Menten kinetic equation for sucrose transport which 
was directly connected to the intracellular Sue pool. The model 
structure is shown in Figure 1. Based on the physical and kinetic 
evidence for an association of SPS and sucrose phosphate phos- 
phatase (SPP; Echeverria et al., 1997), it was assumed that free S6P 
was not occurring, and the interconversion of S6P to Sue was sim- 
ulated in one step with the reaction catalyzed by SPS (indicated 
by dashed line in Figure 1). Rates of metabolite interconversion 
were modeled by Michaelis-Menten kinetics except for the reac- 
tions catalyzed by PGI, PGM, UGPase, TPP and sucrose synthase 
(SuSy). These reaction rates were assumed to be not limiting and 
to follow a mass-action kinetic, i.e., to be directly proportional 
to the substrate concentration. The reaction rate of sucrose cleav- 
age, catalyzed by the invertase (Inv) enzyme, was modeled as an 
irreversible Michaelis-Menten kinetic incorporating mixed inhi- 
bition by the reaction products Glc (non-competitive inhibitor) 
and Fru (competitive inhibitor; Nagele et al., 2010). Reaction rates 
of hexose phosphorylation, catalyzed by the enzyme hexokinase 
(Hxk) with a fructokinase and glucokinase activity, were modeled 
as Michaelis-Menten enzyme kinetics non-competitively inhib- 
ited by F6P (fructokinase) and G6P (glucokinase; Claeyssen and 



Frontiers in Plant Science | Plant Systems Biology 



July 2014 I Volume 5 | Article 365 | 8 



Nagele and Weckwerth 



Linking sugar metabolism to energy homeostasis 



Rivoal, 2007). All other reactions were modeled by irreversible 
non-inhibited Michaelis-Menten kinetics. 

INTEGRATION OF SnRKI ACTIVITY IN SPS KINETICS 

To kinetically integrate the phosphorylation and inactivation of 
SPS by SnRKI, the maximum enzyme velocity of SPS, Vmax,SPS) 
was described as the product of the rate constant /csps and the 
absolute amount of dephosphorylated and, hence, activated SPS 
peptide, £sPS,act (Eq. 2): 

Vmax,SPS = ^SPS ' £sPS,act (2) 

The rate constant ksps and the total enzyme level of SPS 
-Esps,totab i-e., the sum of £sPS,act and £sPS,inact> were estimated from 
previous publications (Osuna etal., 2007; Lehmann etal., 2008). 
We estimated Vmax,SPS = 422 nmol gFW~^ min~^ = 25.32 |xmol 
gFW-^ h-^ (Osuna etal, 2007), and £sPS,total = 8.5 pmol 
gFW-^ = 8.5 X 10"^ ixmolgFW"^ (Lehmann etal, 2008). Based 
on these data we estimated the rate constant ksps (Eq. 3): 

7 1^max,SPS ^ ^ /o\ 

kvs = = 2.98 X 10^ . - (3) 

^SPS,total rl 

Then, v^^^cSPS was integrated in a bi- substrate kinetic to 
describe the reaction rate of SPS (Eq. 4): 

[F6P] [UDPG] 

VSPS Vn,ax,SPS ' ^^^^^^ ^ ^^^^^ ' ^^^^^^^ ^ [UDPG] 

Dynamics in levels of activated enzyme £sps,act and of inac- 
tivated enzyme £sPS,inact were described by ordinary differential 
equations (ODEs) depending on the phosphorylation by SnRKI 
and the dephosphorylation by a protein phosphatase PP (Eq. 5): 

d[£^SPS,act]/d^ = - VSnRK 
d[£^SPS,inact]/df = VSnRK - Vpp 

The rate of SPS activation by dephosphorylation was modeled 
via an irreversible Michaelis-Menten kinetic in which £sps. inact was 
the substrate. The rate of SPS inactivation by phosphorylation was 
modeled by a SnRKI kinetic equation in which £sPS,act was the 
substrate. The derivation of this kinetic equation is described in 
the following section. 

DERIVING A KINETIC EQUATION FOR SnRKI ACTIVITY 

The maximum phosphorylation rate catalyzed by SnRKI was 
estimated based on data from Zhang etal. (2009). The authors 
describe seedling SnRKI activity without the addition of inhibitors 
to be about 2.25 nmol phosphate incorportated in AMARA pep- 
tide min~^ mgProtein"^ (Zhang etal., 2009) which we estimated 
to result in 1.35 |xmol h~^ gFW~^ assuming a protein content 
of about 10 mg gFW~^ (~1% w/w). As reported by Nunes etal. 
(2013b) replacing AMARA peptides by SPS did not significantly 
affect kinetic parameters. Hence, we assumed that these data on 
AMARA peptide phosphorylation are representative for our mod- 
eling approach which focused the SPS phosphorylation. Applying 
data on inhibitory constants for T6P (Ki,T6P = 0.005 mM), GIP 



(Ki,Gip = 0.48 mM) and G6P (Ki,G6P> 1 mM; Nunes etal, 
2013b) as well as the finding that experimental data were fit- 
ted best by a partial non-competitive mixed-type inhibition for 
T6P (Zhang etal, 2009), we then fitted the affinity, i.e., the dis- 
sociation constant KM,SnRKi of SnRKI to SPS peptide using a 
downhill simplex method in multidimensions as it is implemented 
in the Matlab® -based Systems Biology Toolbox 2 (Schmidt and 
Jirstrand, 2006). Simultaneously and applying the same method- 
ology, the ratio of rate constants for partial inhibition by T6P, 
^partial was estimated according to the inhibition kinetics pre- 
sented by Nunes etal.(2013b). Best solutions were obtained for 

J^M,SnrKl = 10"^ mM and /Cpartial = 0-29. 

ODE PROGRAMMING, PARAMETER ESTIMATION AND SENSITIVITY 
ANALYSIS 

The metabolic network shown in Figure 1 was mathemati- 
cally represented in a system of ODEs. The ODEs described 
the time dependent change in metabolite and enzyme levels 
as a function of enzymatic interconversion. Model parameters 
were estimated to simulate changes in metabolite levels during 
a 3 h supply of starved Arabidopsis seedlings with externally 
supplied sucrose with a concentration of 15 mM (Yadav etal., 
2014). Kinetic parameters were numerically estimated by a down- 
hill simplex method in multidimensions which is implemented 
in the Matlab® -based Systems Biology Toolbox 2 (Schmidt and 
Jirstrand, 2006). Upper and lower bounds for kinetic param- 
eters were defined according to literature data (Osuna etal, 
2007; Nagele etal, 2010, 2012). For reactions of PGI, PGM, and 
UGPase no explicit enzyme parameters were available. Instead, 
the reaction catalyzed by cFBPase, which was defined as the input 
reaction of the system, was considered to be the rate limiting 
step (Strand etal, 2000), i.e., resulting rates of PGI-, PGM,- 
and UGPase -driven reactions had to be lower than/equal to the 
maximum enzyme activity of cFBPase which was analyzed by 
Osuna etal. (2007). The resulting ODE model is provided with 
parameters to simulate dynamics of carbohydrate metabolism 
induced by external sucrose supply (Supplementary Information 
1: Model_dynamic.txt) as well as to simulate a metabolic steady- 
state after 3 h of sucrose supply (Supplementary Information 2: 
Model_steady-state.txt) . The model can be simulated and analyzed 
within the Matlab® -based Systems Biology Toolbox 2 (Schmidt 
and Jirstrand, 2006). 

Sensitivity analysis was performed using the "SBsensdatastat" 
and the "SBsensstat" function implemented in the Systems Biology 
Toolbox 2 (Schmidt and Jirstrand, 2006). During this sensitivity 
analysis, model parameters are perturbed with a defined intensity 
(here: 1% of the parameter value). Afterward, simulation results of 
the perturbed system are compared to the non-perturbed system. 
In this way it is possible to estimate the impact of small param- 
eter perturbations to model states, i.e., metabolite and protein 
concentrations and reaction rates. 
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